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Abstract. The size estimates approach for Electrical Impedance Tomography 
(EIT) allows for estimating the size (area or volume) of an unknown inclusion 
in an electrical conductor by means of one pair of boundary measurements of 
voltage and current. In this paper we show by numerical simulations how to 
obtain such bounds for practical application of the method. The computations 
are carried out both in a 2-D and a 3-D setting. 



1. Introduction 

EIT is aimed at imaging the internal conductivity of a body from current and 
voltage measurements taken at the boundary. It is well known, [I], P], that, even in 
the ideal situation in which all possible boundary measurements are available, the 
correspondence boundary data — > conductivity is highly (exponentially) unstable. 
As a consequence it is evident that, in practice, it is impossible to distinguish high 
resolution features of the interior from limited and noisy boundary data, [3j . 

Motivated by applications, a line of investigation pursued by many authors, [4], 
[5], [6], [7], 0, [9], [10], [11], has been the one of hmiting the analysis to cases in 
which one seeks an unknown interior inclusion embedded in an otherwise known 
(may be even homogeneous) conductor, and whose conductivity is assumed to differ 
from the background. 

Even in this restricted case, and even when full boundary data are available, the 
instability remains of exponential type [12j . 

It is therefore reasonable to further restrict the goal and attempt to evaluate 
some parameters expressing the size (area, volume) of the inclusion, disregarding 
its precise location and shape, having at our disposal one pair of boundary mea- 
surements of voltage and current. This approach, which can be traced back to [4], 
has been well developed theoretically, [13], [M], [15], [16], see also [17] and [H] for 
the analogous treatment in the linear elasticity framework. In order to describe 
such type of results we need first to introduce some notation. 

We denote by a bounded domain in M", n — 2,3, representing an electrical 
conductor. The boundary d^l of fl is assumed of Lipschitz class, with constants 
To, Mq, that is the boundary can be locally represented as a graph of a Lipschitz 
continuous function with Lipschitz constant Mq in some ball of radius rg . When no 
inclusion is present in the conductor we assume that it is homogeneous and we pose 
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its conductivity (7{x) = 1. When the conductor contains an unknown inclusion D 
of different conductivity, say fc > 0, fc ^ 1 the overaU conductivity in the conductor 
wiU be given by cr{x) = 1 + (fc — 1)xd{x). Here and in what foUows it is assumed 
that D is strictly contained in ^l. More precisely, for a given do > 0, 

dist{D,dn) >do. (1.1) 

Let (f G H^^{dVL), Jg^f = 0, be an applied current density on dft. The induced 
electrostatic potential u £ H^{n) is the solution of the Neumann problem 

div ((1 + (fc - 1)xd) Vu) = 0, in n, 

(1.2) 

Vu ■ 1^ = if, on dfl, 

where i' denotes the outer unit normal to dft. 

When D is the empty set, that is when the inclusion is absent, the reference 
electrostatic potential uq £ iJ^(ri) satisfies the Neumann problem 

Auq — 0, in ri, 

(1.3) 

Vuq ' ^ = on d^- 
In both cases (jl.2p and (jl.3p . the solutions u and uq are determined up to an 
additive constant. 

Let us denote by VF, Wq the powers required to maintain the current density ip 
on dfl when D is present or it is absent, respectively. Namely 

f uifi^ f (l + (fc-l)xi3)|Vup, (1.4) 

Wo^ f uoifi^ f iVuop. (1.5) 
Jdfi Jn 

The size estimate approach developed in [TJ], [T3], [TS], [T^, tells us that the 
measure \D\ oi D can be bounded from above and below in terms of the quantity 

^wd^ which we call the normalized power gap. More precisely the following 

bounds hold, see fW, Theorem 2.3]. 

Theorem 1.1. Let D be any measurable subset of satisfying (jl.ip . Under the 
above assumptions, if k > 1 we have 

Ct%r^<\D\< -—) C+ . (1.6) 



k-1 ' Wo - ' ' - \k-lj ^ V 
If conversely, k < 1, then we have 

':Or^^im^(^'c,(^^\ (1.7, 



1 - fc ^ Wo - ' ' - \l-kj ^ V Wo 

where , Ci only depend on do, ro. Mo, whereas p > 1, C2 , C2 only depend 
on the same quantities and, in addition, on the frequency of (p 

FW] = ^^iiim. (1.8) 

When it is a priori known that the inclusion D is not too small (if it is at all 
present), a situation which often occurs in practical applications, stronger bounds 
apply. 
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Theorem 1.2. Under the above hypotheses, let us assume, in addition, that 

\D\>mo, (1.9) 
for a given positive constant mg. If k > 1 we have 

fc - 1 1 Wo - - A: - 1 2 Wo ^ ' 

If, conversely, k < I, then we have 

^ .cr2^ < 1^1 < ^ci^. (1.11) 



l~k ' Wo - ' ' - 1 - k ^ Wo 

where , Cj" only depend on do, \n\, ro, Mo, whereas C2 , C2 only depend on 
the same quantities and, in addition, on nio and F[Lp\. 



Theorem 1 1.2 1 can be easily deduced from Theorem 1 1.1 1 by the arguments sketched 
in [191 Appendix]. 

One of the goals of the present paper is to test the applicability of such bounds 
by numerical simulations with the following purposes: 

i) provide practical evaluations of the constants , appearing in the above 
inequalities (fTTU)) . pTTT]) .- 

ii) when, due to special geometric configurations, it is possible to compute the- 
oretically such constants, compare such theoretical values with those obtained by 
simulations; 

Hi) show that such upper and lower bounds deteriorate as the frequency F[ip\ 
increases. 

The other goal of this paper is to perform similar kinds of numerical simulations 
when the so-called complete model of EIT is adopted. We recall that this model 
is aimed at an accurate description of the boundary measurements suitable for 
medical applications, and was introduced in 20J and subsequently developed in 
PT] and . In this model, the metal electrodes behave as perfect conductors and 
provide a low-resistance path for current. An electrochemical effect at the contact 
between the electrodes and the body results in a thin, highly resistive, layer. The 
impedance of this layer is characterized by a positive quantity zi on each electrode 
e;, I = 1, L, which is called surface impedance. 

Denoting by /; the current applied to each e/, the resulting boundary condition 
on each electrode e; becomes 

u + ziVu-u = U\ one;, (1.12) 

where the unknown constant [/' is the voltage which can be measured at the elec- 
trode e;. 

We assume, as before, that the reference conductor has conductivity a = 1 
and that an unknown inclusion D of conductivity a = k, with fc > and fc 7^ 1, is 
strictly contained in fi. Therefore, the electrostatic potential u inside the conductor 
is determined, up to an additive constant, as the solution to the following problem 

div((l + (fc-l)xD)VM) = 0, inQ, 

u + ziVu-v ^U'', onei,l<l<L, ciT^A 

Vu-v = 0, ondn\utiei, 

J Vu -1^ = 11, 1<1<L, 
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where the so-cahed current pattern / = (/i, /l) is subject to the conservation of 
charge condition ~ 0, and the unknown constants U' are the components 

of the so-caUed voltage pattern U = {U^, U^). 

When the inclusion is absent, the electrostatic potential uq induced by the same 
current pattern / is determined, up to an additive constant, as the solution of the 
following problem 

Auq = 0, in fi, 

uo + ziVuo ■ = Ul^, one/, l</<i, 
' Vuo-i'^O, ondn\utiei, ^ ' 

^ IVu^-v^h, 1<1<L, 

where, as before, the Uq are unknown constants in the direct problem p.l4p . 

We shall assume that the sets ei,...,ei, representing the electrodes, are open, 
pairwise disjoint, connected subsets of dfl and, in addition, 

dist(ei, Cfc) > ^1 > for every /, fc, I ^ k. (1-15) 

The surface impedance on e;, / = 1, L, is assumed to be real valued and to 
satisfy the following bounds 

< m < Zi < M, for every I = 1, L. (1.16) 

In this formulation, the powers W and Wq become 

L 

VF = ^/,{/\ (1.17) 

1=1 

L 
i=l 

Size estimates like those of Theorems 1 1 . 1 [ fTT2] were obtained for the complete model 
in [23] . In particular we have 

Theorem 1.3. Let D be any measurable subset of fl satisfying (jl.ip and let W, 
Wq be given by (jl.l7p . (ll.lSp . Then, inequalities (|1.6p . (jl.7p hold for fc > 1 and 
fc < 1, respectively, where the constants , Ci only depend on do, rg, Mq, 
and C2 , C2 and p > 1 only depend on the same quantities and, in addition, on 5i, 
M and m. 

Also in this case, the size estimates of \D\ can be improved to the form (jl.lOp . 
(|l.lip when condition (|1.9p is satisfied. 

In Section [5] we consider the standard EIT setting. We start by describing the 
finite element setup used in our numerical simulations in Section 12.11 Next (as a 
warmup) we consider a two-dimensional model in Section 12.21 

In Section [2T3I we consider the three-dimensional case and we discuss all items i), 
ii), iii) introduced above. In particular we observe that, comparing the results as 
the frequency F[(p] increases, we have quite rapidly a serious deterioration of the 
bounds. This poses a severe warning on the limitations that have to be taken into 
account in the choice of the boundary current profile if. 

Section [3] is devoted to simulations with the complete EIT model. In this case 
it is reasonable analyze the case when only two electrodes, one positive and one 
negative, are attached to the surface of the conductor. In this case, the frequency 
function is not available from the data since we are not prescribing the boundary 
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current Vm • v\dn but only the current pattern, which is a 2-electrode configuration, 
is just the pair (1, —1). In place of the frequency function, the parameters that may 
influence the constants in the volume bounds are: the width of the electrodes and 
the distance between them. We perform various experiments to test such variability. 



2. Numerical simulations for the EIT model 

2.1. Numerical model. The numerical model is based on the discretization of the 
energy functional J : ffi(17,K") M 

J{u) = \ ! {l + {k- l)xD)yu -Vu- f ^u, (2.1) 

associated to the variational formulation of problem (jl.2[) . The energy functional 
(12. ip has been discretized by using the High Continuity (HC) technique already 
presented in [21] and [5S] in the context of linear elasticity. Accordingly, for 2-D 
problems the electric potential on the e-th finite element can be represented as 

3 

U,= Y,'P^{^M^2)U^,, (2.2) 

whereas for the 3-D case it assumes the form 

3 

Ue = ^ M^l)4>j{^2)(t)l{£.3)Utjl, (2.3) 
i.j,l=l 

where the coordinates £,r, ^ = 1, ■■■,n, span the unitary element domain [— ^, ^]", 
n = 2, 3, and Uij, u^i are the HC parameters involved in the field interpolation on 
the generic element. The shape functions 0i(^r) are defined as 





1 

^8 


2^ 








3 

~4 






(2.4) 


V 


1 

^8 









The 1-D case illustrated in Figure [T] shows the meaning of the HC parameters. 
They allow to define the slopes of the interpolated function at the end points of the 
element. On the same figure one can see also the positions of the HC nodes and 
the shape functions (|2.4p . 

Figure [5] shows a typical structured mesh on a rectangular domain and the nodes 
used for the approximation of the potential field in the 2-D case. For elements with 
a side lying on the boundary, in order to easily impose the Neumann boundary 
conditions, special shape functions are used. In practice, the external HC nodes 
are translated onto the boundary dO. and the related HC parameters have the 
meaning of function values (see again Figure [5]). In this case the shape functions 
relative to a left boundary (^^ — —\) and a right boundary (^^ = \) of the finite 
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Figure 1. HC interpolation in the 1-D case: nodes, parameters 
and shape functions. 




Figure 2. HC mesh in the 2-D case: nodes for boundary and 
inner elements. 



element are 
left : < 02 (Cr 

03 (Cr 



-e- 



right : < 



03 (^r 



(2.5) 



Further details about the HC interpolation can be found in 24! and |25]. This 
interpolation technique, which can be considered as a particular case of the Bezier 
interpolation, has the main advantage of reproducing potential fields of smooth- 
ness with a computational cost equivalent to a interpolation. 

By (I2.2p or (|2.3[) , the potential field u on each element e takes the compact form 

Ue = NeWe. (2.6) 

The one-row matrix Ng collects the shape functions of the HC interpolation, 
whereas the components of the vector Wg are the nodal parameters of the un- 
derlying element. With this notation, the gradient of the potential field is given 

by 



(2.7) 
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We remark that the dimensions of the matrices Ne, VNe and vector We are 
1 X 9, 2 X 9 and 9 X 1 for the 2-D case and 1 x 27, 3 x 27 and 27 x 1 for the 3-D 
case. 

By (|2.6p and (|2.7p . the discrete form of (|2.ip becomes 

J(We) = X] Q (1 + (fc - l)Xl?)(VNeWe) • (VN^We) - ^^NeW^^ , (2.8) 

or, in a compact form, 

J(W,) = J2 (wfKeWe - pe) . (2.9) 

e 

The latter equation provides the definition of the matrix and vector associated to 
e-th element 

' Ke = /" (1 + (fc - l)xD)(VNe)^VNe, 

J?= (2.10) 

Jdn^ 

which can be used to assemble, by using standard techniques, the system of equa- 
tions to solve. 

2.2. Tvifo— dimensional case. Numerical analysis has been performed on a square 
conductor Q of side I under the two current density fields ip illustrated in Figure [3l 
The domain has been discretized with a mesh of 21 x 21 HC finite elements and 
for both Test Ti and Test T2 of Figure [3] we have considered an inclusion D with 
conductivity fc = 0.1orfc = 10. 



(p=0 



^=1/4 




(p=0 
(a) 



Figure 3. Square conductor considered in 2-D numerical simu- 
lations for the EIT model and applied current density fields: Test 
Ti (a). Test T2 (b). 



A first series of experiments has been carried out by considering all possible 
square inclusions with side ranging from 1 to 5 elements, that is the size of inclusion 
has been kept lower than 6% of the total size of the conductor. The results are 
collected in Figures |4] and [5] for different values of the minimum distance do between 
the inclusion D and the boundary of f2. 

From Figures m^a) and[5fa), which refer to the case k = 0.1, one can note that 
the upper bound of \D\ is rather insensitive to the choice of do, whereas the lower 
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0.00 0.05 0.10 0.15 0.20 0.25 0.00 0.05 0.10 0.15 0.20 

(W-Wnj/Wo {W„-W)IWo 



(a) 



(b) 



Figure 4. Influence of do for square inclusions in Test Ti of Figure 
Elja) (21 X 21 FE mesh): k = 0.1 (a), fc = 10 (b). 




0.00 0.05 0.10 0.15 0.20 0.25 0.00 0.05 0.10 0.15 0.20 0.25 

(»-H'„)/H'„ (H'o-HO/W,, 



(a) 



(b) 



Figure 5 . Influence of do for square inclusions in Test T2 of Figure 
ini;b) (21 X 21 FE mesh): k = 0.1 (a), fc = 10 (b). 



bound in (|l.lip improves as do increases. The converse situation occurs when the 
inclusion is made by material of higher conductivity, see Figures lUJb) and Eljb) . 

As a second class of experiments, we have considered inclusions of general shape 
on a FE mesh of 15 x 15 HC elements. More precisely, each inclusion is the union 
of elements having at least a common side and being at least do = 2 elements far 
from the boundary dO.. Results are collected in Figures [5] and [71 

The straight lines drawn in Figures [4] and [6] correspond to the theoretical size 
estimates for test Ti of Figure [3l^a). For both cases fc = 0.1 and fc = 10 we have 

l \W-W^\ ^ \D\ ^ 10 |Vt^-T4^o| , . 

9 Wo - \n\ - % Wq ' ^ ' ' 
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The comparison with the region of the plane - — wiT^^J covered by the corre- 
sponding numerical experiments confirms, as already remarked in [19| in the context 
of linear elasticity, that practical applications of the size estimates approach lead to 
less pessimistic results with respect to those obtained via the theoretical analysis. 




(a) (b) 

Figure 6. Numerical size estimates for inclusions of general shape 
generated from a generic element inside ft for test Ti of Figure [3ja) 
(21 X 21 FE mesh, do = 2): fc = 0.1 (a), fc = 10 (b). 




Figure 7. Numerical size estimates for inclusions of general shape 
generated from a generic element inside D, for test T2 of Figure[31Jb) 
(21 X 21 FE mesh, do = 2): k = 0.1 (a), fc = 10 (b). 
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2.3. Three— dimensional case. The first part of this subsection is devoted to the 
extension to the 3-D case of the numerical simulations given in 12.21 In the second 
part, we shall investigate on the effect of the oscillation character of the Neumann 
data on the upper bound of size inclusion. 

Similarly to the 2-D first series of numerical simulations has been per- 

formed on an electrical conductor of cubic shape, of side I, with the two current 
density fields illustrated in Figure El In both mesh of 20 x 20 x 20 finite 

elements has been considered when performing simulations in presence of cubic in- 
clusions. The results are illustrated in Figures and [TUl Figure contains also the 
straight lines corresponding to the theoretical size estimates for test Ti of Figure 
El that is 

l \W-Wo\ ^ \D\ ^ lO lTy-T^ol 
9 T^o " " 9 1^0 " 






^=1/8 






93=1/8 


^=1/2 




V3=l/8 




^=1/8 






(b) 





(a) 



Figure 8. Cubic conductor considered in 3-D numerical simula- 
tions for the FIT model and applied current density fields: Test Ti 
(a) and Test Ta (b). 



In order to deal with inclusions of general shape, however, the numerical exper- 
iments require some restrictions to reduce the computer time. A rough estimate of 
the computational cost can be obtained noting that the numerical effort is essen- 
tially due to the decomposition of the matrix associated to the linear system (12. 9p 
and to the computation of its solution. Denoting by m the number of the equations 
and by b the half bandwidth of the matrix, the decomposition requires m{b — 1) 
multiplications and m&(6 — 1) additions, whereas the computation of the solution 
involves mb multiplications. 

Therefore, for each given inclusion in a 20 x 20 x 20 FE mesh, a linear system of 
10648 (6 — 1015) equations has to be solved, requiring a computer time of approx- 
imately 86 s working on an Opteron 2.4 GHz computer. Since the number of all 
possible inclusions formed by Ui elements on a mesh of rie x rie x rie is „ : 
the way to calculate all the possible case is practically impossible. Indeed by con- 
sidering that the 20 x 20 x 20 is formed by 8000 elements and that, if the ratio 
|D|/|ri| is less than 6% that is 480 elements, the number of cases to analyze is 
69.1183 X lO''^^ 
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Figure 10. Influence of do for cubic inclusions in test T2 of Figure 
E (20 X 20 X 20 FE mesh): k = 0.1 (a), A: = 10 (b). 



In order to reduce the computer time significantly we have considered a 7 x 7 x 7 
mesh generating a system of 729 equations. Despite of this, the number of possible 
cases to consider still remains very high; for instance, for inclusions formed by 5 
elements, one should solve about 3.8 x 10^° linear systems. Therefore, we decided 
to restrict our analysis to inclusions satisfying the following additional hypotheses: 

i) the inclusion is the union of elements having at least one common face and it 
is formed by starting from a generic element inside an octant of the cube (this 
last assumption is not really restrictive due to the symmetries of the problem) ; 

ii) do = 1. 
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For inclusions formed by 1, 7 elements, we have considered all possible inclusions 
satisfying the limitations i) and ii), whereas for inclusions formed by 8, 17 ele- 
ments we have considered a random sample because of the high computational cost. 
For these cases, the ratio between the sample dimension and that of all the data 
approximately spans between 20% for inclusions formed by 8 elements and 0.01% 
for inclusions formed by 17 elements. The results are presented in Figures fTT] and [T2l 
for Test Ti and Test T2, respectively. In Figure [TTJ the straight lines corresponding 
to the theoretical bounds (|2.12p for Test Ti are also drawn. As already remarked in 
the treatment of the 2-D case, the theoretical analysis leads to rather pessimistic 
results with respect to those obtained by the numerical simulations, especially when 
the inclusion is softer than the surrounding material. 




(a) (b) 

Figure 11. Numerical size estimates for inclusions of general 
shape generated from a generic element belonging to an eight of 
the cube for test Ti of Figure [ija) (7 x 7 x 7 FE mesh, do = 1): 
k = 0.1 (a), k^lO (b). 

The Neumann data considered in the above experiments give raise to potential 
fields inside the conductor with nonvanishing gradient. In the general case, when 
the gradient of the solution may vanish, we expect, accordingly to Theorems 11.11 [L^ 
that the upper bounds deteriorate as the frequency F[ip\ given by (jl.8p increases. 
Since F[ip] is a ratio which measures the frequency of oscillation of ip, we are 
interested to investigate on the effectiveness of size estimates approach for oscillating 
Neumann data. 

In particular, the numerical simulations have been carried out for the cubic 
electrical conductor considered before and choosing the following Neumann data: 

(fi — - cos on z = 0, 1 

ip^ cos^ onz^l, \ torn = 0,1, 2. (2.13) 

(fi = elsewhere on 9ri, J 

Case n = has been already discussed at the beginning of this paragraph and 
corresponds to the simple case in which the gradient of the unperturbed solution 
uo does not vanish in fl. 
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Figure 12. Numerical size estimates for inclusions of general 
shape generated from a generic element belonging to an eight of 
the cube for test T2 of Figure [H; a) (7 x 7 x 7 FE mesh, do = 1): 
k = 0.1 (a), /c = 10 (b). 



The two other cases are examples of Neumann data with higher frequency F[lp\ 
More precisely, the corresponding solutions mq have critical lines of equation 



+ i 



I (I 



= 0,, 



1. 



The mesh employed is made by 20 x 20 x 20 HC finite elements. The analysis 
has been focussed on cubic inclusions having volume up to 6% of the total volume 
of the specimen and conductivity fc = 0.1 and fc — 10. The numerical results in case 
n = 1 and n = 2 are presented in Figures [13] and 1141 respectively. The numerical 
results show that the lower bound in size estimates (jl.lOp . (|1.11|1 improves as do 
increases, whereas the upper bound of jZ?! is rather insensitive to the choice of do. 

Theoretical estimates for cases n = 1 and rt = 2 of (|2.13p are given by 

for fc > 1 : 

-W \D\ Ik 
1) W^-joy-afc^' ^2.14) 



tanh^f^ 



mr{k — 
for fc < 1 : 



tanh '-f Wo 
riTT 



W 



Wo 



n7r(l — k) 



tanh 



mrW-Wo 



Wo 



< 



\£1 



< 



1 1 tanh ifW-Wo 



Cn 1 k 



Wo 



where 



10 



rnr cosh 



f sinh — 
V 20 



sm — 
20/ 



> 7^=1,2. 



The theoretical estimates are indicated in Figures [T3] and [TH The slope of the 
straight line corresponding to the upper bound is so high that it practically coincides 
with the vertical axis, at least for the portion of graph near the origin considered 
in this study. The theoretical lower bound gives, for a fixed power gap, values 
significantly less than those obtained in the numerical experiments. 
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(a) (b) 

Figure 13. Cubic electrical conductor with Neumann data as in 
case n = 1 of (|2.13p : lower and upper bound of the power gap for 
different values of do (fc = 0.1 (a) and fc = 10 (b)) on a 20 x 20 x 20 
mesh. 




0.00 0.05 0.10 0.15 0.20 0.25 0.00 0.05 0.10 0.15 0.20 0.25 

{W-Wo)IWa {Wa-W)IWo 



(a) (b) 

Figure 14. Cubic electrical conductor with Neumann data as in 
case n = 2 of (j2.13p : lower and upper bound of the power gap for 
different values of do (fc = 0.1 (a) and fc = 10 (b)) on a 20 x 20 x 20 
mesh. 



3. Numerical simulations for the complete EIT model 



3.1. Numerical model. In this case, by using the same notation introduced in 
Section [21 the energy functional J : H^iVl) x — > R related to the variational 
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formulation of problem (|1.13l) is given by 

J{u, U') = l fil + {k- 1)xd)Vu .Vu+IY.- [ - f^')' - E (3.1) 

Using HC interpolation for the potential field u and with the notation introduced 
in Section [21 the discrete energy functional becomes 



J(We,i7') = i^wfKeWe + 



(3.2) 



1=1 e i=l 

having used the compact notation 

Ke = / (1 + (fc - l)XD)(VN,)^VNe, 

Ku = I NjN, 



(3.3) 



(ao,)e 



(3.4) 



N 



(ao,)e 



We remark that the second sum in the right hand side of (|3.2p and (|3.3[) . that on 
e, is extended only to the elements under the electrodes. 

Collecting the unknown parameters representing the potential field in w, those 
of the electrodes in U and the current pattern in I, by a standard method of 
assembling we obtain the following linear system 





w 




o" 




U 




I 



(3.5) 



which can be efficiently solved taking advantage of the particular structure of coef- 
ficient matrix. 

3.2. Results for 3— D cases. The analysis has been restricted to the case of two 
electrodes located on the boundary of a cubic electrical conductor of side I, see 
Figure [Tsl The specimen has been discretized by a mesh ofl7xl7xf7 cubic 
HC finite elements and the numerical experiments have been carried out on cubic 
inclusions only, with volume up to 6% of the total volume and conductivity value 
k = Q.l or k = fO. The surface impedance takes a constant value such that 
^ = ^ = 0.2 on both electrodes, according to properties of human skin reported 
in literature, see, for instance, [55]. 

In test Ti of Figure [151 the electrodes cover completely two opposite faces of the 
specimen, whereas in Test T2 one electrode coincides with a face of dft and the 
other is a square, formed by one or nine surface finite elements, and it is located in 
central position of the opposite face. Finally, in Test T3, two electrodes are placed 
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on the same face of the conductor in a symmetric way respect to middle hnes of 
the face. The electrodes are separated by three finite elements and their dimensions 
are equal to the element size. 
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(a) (b) (c) 



Figure 15. Cubic conductor considered in 3-D numerical simu- 
lations for the physical EIT model and location of the electrodes: 
test Ti (a), test (b) and test (c). 



The numerical results for Test Ti are presented in Figure [11] for fc 
fc = 10, respectively, and for varying values of dg- For both cases k 
fc = 10, the theoretical size estimates are given by 
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and, again, they lead to a rather pessimistic evaluation of the upper and lower 
bounds. 
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Figure 16. Influence of do for cubic inclusions in Test Ti of Figure 
EJa) (17 X 17 X 17 FE mesh, C = 0.2 ): 0.1 (a). A: = 10 (b). 



Concerning Test T2, Figure fTTl shows the results when the small electrode coin- 
cides with one surface finite element, whereas Figure [18] refers to the case of a 3 x 3 
finite elements electrode. One can notice that in all the four cases, the upper bound 
is not really influenced by the value of dg. Moreover, the inaccuracy in determining 
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the lower bound of the angular sector, is probably due to the fact that the present 
analysis is restricted to the special class of cubic inclusions. 

A comparison between Figure [T7] and Figure [TS] suggests that better upper 
bounds can be obtained by enlarging the size of the small electrode. Moreover, 
from Figures [17] and [18] it appears clearly that the lower bound significantly im- 
proves as the distance do between the inclusion D and the boundary of f2 increases. 
This property has been further investigated by increasing only the distance dos 
of the inclusion D from the face of the conductor containing the small electrode. 
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Figure [19] shows the resuhs of simulations in the case of a single finite element 
electrode and a comparison with Figure [17] suggests that the improvement of the 
lower bound is mainly due to the greater distance from the electrode. 




0.00 0.10 0.20 0.30 0.40 0.50 0.00 0.10 0.20 0.30 0.40 0.50 

{W-Wo)IWo {Wa-W)IW„ 



(a) (b) 

Figure 19. Influence of doa for cubic inclusions in Test T2 of Fig- 
ure [T5i;b) (17 X 17 X 17 FE mesh, C = 0.2, 1 x 1 FE electrode): 
k = 0.1 (a), /c = 10 (b). 




Figure 20. Influence of for cubic inclusions in Test T3 of Figure 
P^c) (17 X 17 X 17 FE mesh, C = 0.2, 1 x 1 FE electrode): k = 0.1 
(a), fc = 10 (b). 



Finally, the results of the numerical simulations for Test T3 are presented in 
Figure [20] In this case, the lower bound improves as the distance do between 
the inclusion D and the boundary of 17 increases, whereas the upper bound is 
indistinguishable from the vertical axis. 
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4. Conclusions 

We have tested by numerical simulations the approach of size estimates for EIT. 
We could perform experiments in the 2-D setting with a large varieties of shapes of 
inclusions and we found quite satisfactory bounds, which in some cases are markedly 
better than those derived theoretically. 

In the 3-D case, we had to limit the variety of shapes of the test inclusions since 
the growth of their degree of freedom conflicts with the limitations on computer 
time. We showed that good volume bounds hold when the boundary data is well- 
behaved in terms of its frequency, whereas they rapidly deteriorate as the frequency 
increases. 

For the complete EIT model we have also made tests in a 3-D setting and 
compared the bounds in terms of the size of the electrodes, their relative distance 
and their a-priori assumed distance from the inclusion D. We have shown that we 
obtain good bounds when the electrodes are not too small and when D is sufficiently 
away from them. 
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